The limits of mean-field heterozygosity estimates under spatial extension in simulated plant populations

Computational models of evolutionary processes are increasingly required to incorporate multiple and diverse sources of data. A popular feature to include in population genetics models is spatial extension, which reflects more accurately natural populations than does a mean field approach. However, such models necessarily violate the mean field assumptions of classical population genetics, as do natural populations in the real world. Recently, it has been questioned whether classical approaches are truly applicable to the real world. Individual based models (IBM) are a powerful and versatile approach to achieve integration in models. In this study an IBM was used to examine how populations of plants deviate from classical expectations under spatial extension. Populations of plants that used three different mating strategies were placed in a range of arena sizes giving crowded to sparse occupation densities. Using a measure of population density, the pollen communication distance (Pcd), the deviation exhibited by outbreeding populations differed from classical mean field expectations by less than 5% when Pcd was less than 1, and over this threshold value the deviation significantly increased. Populations with an intermediate mating strategy did not have such a threshold and deviated directly with increasing isolation between individuals. Populations with a selfing strategy were influenced more by the mating strategy than by increased isolation. In all cases pollen dispersal was more influential than seed dispersal. The IBM model showed that mean field calculations can be reasonably applied to natural outbreeding plant populations that occur at a density in which individuals are less than the average pollen dispersal distance from their neighbors

Charge environments around phosphorylation sites in proteins

Background: Phosphorylation is a central feature in many biological processes. Structural analyses have identified the importance of charge-charge interactions, for example mediating phosphorylation-driven allosteric change and protein binding to phosphopeptides. Here, we examine computationally the prevalence of charge stabilisation around phosphorylated sites in the structural database, through comparison with locations that are not phosphorylated in the same structures. Results: A significant fraction of phosphorylated sites appear to be electrostatically stabilised, largely through interaction with sidechains. Some examples of stabilisation across a subunit interface are evident from calculations with biological units. When considering the immediately surrounding environment, in many cases favourable interactions are only apparent after conformational change that accompanies phosphorylation. A simple calculation of potential interactions at longer-range, applied to non-phosphorylated structures, recovers the separation exhibited by phosphorylated structures. In a study of sites in the Phospho.ELM dataset, for which structural annotation is provided by non-phosphorylated proteins, there is little separation of the known phospho-acceptor sites relative to background, even using the wider interaction radius. However, there are differences in the distributions of patch polarity for acceptor and background sites in the Phospho.ELM dataset. Conclusion: In this study, an easy to implement procedure is developed that could contribute to the identification of phospho-acceptor sites associated with charge-charge interactions and conformational change. Since the method gives information about potential anchoring interactions subsequent to phosphorylation, it could be combined with simulations that probe conformational change. Our analysis of the Phospho.ELM dataset also shows evidence for mediation of phosphorylation effects through (i) conformational change associated with making a solvent inaccessible phospho-acceptor site accessible, and (ii) modulation of protein-protein interactions

MCMC-ODPR : primer design optimization using Markov Chain Monte Carlo sampling

Background Next generation sequencing technologies often require numerous primer designs that require good target coverage that can be financially costly. We aimed to develop a system that would implement primer reuse to design degenerate primers that could be designed around SNPs, thus find the fewest necessary primers and the lowest cost whilst maintaining an acceptable coverage and provide a cost effective solution. We have implemented Metropolis-Hastings Markov Chain Monte Carlo for optimizing primer reuse. We call it the Markov Chain Monte Carlo Optimized Degenerate Primer Reuse (MCMC-ODPR) algorithm. Results After repeating the program 1020 times to assess the variance, an average of 17.14% fewer primers were found to be necessary using MCMC-ODPR for an equivalent coverage without implementing primer reuse. The algorithm was able to reuse primers up to five times. We compared MCMC-ODPR with single sequence primer design programs Primer3 and Primer-BLAST and achieved a lower primer cost per amplicon base covered of 0.21 and 0.19 and 0.18 primer nucleotides on three separate gene sequences, respectively. With multiple sequences, MCMC-ODPR achieved a lower cost per base covered of 0.19 than programs BatchPrimer3 and PAMPS, which achieved 0.25 and 0.64 primer nucleotides, respectively. Conclusions MCMC-ODPR is a useful tool for designing primers at various melting temperatures at good target coverage. By combining degeneracy with optimal primer reuse the user may increase coverage of sequences amplified by the designed primers at significantly lower costs. Our analyses showed that overall MCMC-ODPR outperformed the other primer-design programs in our study in terms of cost per covered base

Early jump-out corner detectors

We present two new corner detectors; one works by using dissimilarity along the contour direction to detect curves in the image contour and the other estimates image curvature along the contour direction. These operators are fast, robust to noise, and require no subjective thresholding

Surprisingly low limits of selection in plant domestication

Current debate concerns the pace at which domesticated plants emerged from cultivated wild populations and how many genes were involved. Using an individual-based model, based on the assumptions of Haldane and Maynard Smith, respectively, we estimate that a surprisingly low number of 50–100 loci are the most that could be under selection in a cultivation regime at the selection strengths observed in the archaeological record. This finding is robust to attempts to rescue populations from extinction through selection from high standing genetic variation, gene flow, and the Maynard Smith-based model of threshold selection. Selective sweeps come at a cost, reducing the capacity of plants to adapt to new environments, which may contribute to the explanation of why selective sweeps have not been detected more frequently and why expansion of the agrarian package during the Neolithic was so frequently associated with collapse

The Evolution of Legal Risks Pertaining to Patch Management and Vulnerability Management

This article begins with an overview, in non-technical terms, of the tools generally available and processes implemented for vulnerability management and patch management. Section II identifies some of the evolving security standards that regulators and plaintiffs may rely on to show that companies are legally required to have vulnerability management and patch management. Section III identifies U.S. legal implications of vulnerability management and patch management and factors that a court and regulators may consider

Particle size distributions and the vertical distribution of suspended matter in the upwelling region off Oregon

Various methods of presenting and mathematically describing particle size distribution are explained and evaluated. The hyperbolic distribution is found to be the most practical but the more complex characteristic vector analysis is the most sensitive to changes in the shape of the particle size distributions. Particle size distribution, nutrient concentrations, temperature and other biological and hydrographic data were taken during two cruises off the Oregon coast. The first, in late July, 1973, was during a period of consistent upwelling-favorable winds. The second cruise, in August, 1974, was during a period of intermittently favorable winds. Thus the data presented represent several different upwelling situations. Two distinct vertical structures of suspended particulates and two types of particle size distributions were found, separated by a particle front. On the offshore side of the front, the structure was characterized by dominantly small particles and a subsurface maximum of suspended matter. On the other side of the front, the structure shows a particle maximum at or very near the surface with dominantly large particles. A method for determining onshore-offshore flow patterns from the distribution of particulates was presented. The method was applied to the data from the two cruises. A further experiment was suggested with an emphasis on determining three-dimensional current patterns at the same time as particle distributions. Such data would be used for a three-dimensional numerical model. A numerical model of the vertical structure of two size classes of particles was developed. The results show a close similarity to the observed distributions but overestimate the particle concentration by forty percent. This was attributed to ignoring grazing by zooplankton. Sensitivity analyses showed the size preference was most responsive to the maximum specific growth rates and nutrient half saturation constants. The vertical structure was highly dependent on the eddy diffusivity followed closely by the growth terms

Concentration and δD of molecular hydrogen in boreal forests: Ecosystem-scale systematics of atmospheric H_2

We examined the concentration and δD of atmospheric H2 in a boreal forest in interior Alaska to investigate the systematics of high latitude soil uptake at ecosystem scale. Samples collected during nighttime inversions exhibited vigorous H_2 uptake, with concentration negatively correlated with the concentration of CO_2 (−0.8 to −1.2 ppb H_2 per ppm CO_2) and negatively correlated with δD of H_2. We derived H_2 deposition rates of between 2 to 12 nmol m^(−2) s^(−1). These rates are comparable to those observed in lower latitude ecosystems. We also derive an average fractionation factor, α = D:H_(residual)/D:H_(consumed) = 0.94 ± 0.01 and suggestive evidence that α depends on forest maturity. Our results show that high northern latitude soils are a significant sink of molecular hydrogen indicating that the record of atmospheric H_2 may be sensitive to changes in climate and land use

Elevated glutamatergic compounds in pregenual anterior cingulate in pediatric autism spectrum disorder demonstrated by 1H MRS and 1H MRSI.

Recent research in autism spectrum disorder (ASD) has aroused interest in anterior cingulate cortex and in the neurometabolite glutamate. We report two studies of pregenual anterior cingulate cortex (pACC) in pediatric ASD. First, we acquired in vivo single-voxel proton magnetic resonance spectroscopy ((1)H MRS) in 8 children with ASD and 10 typically developing controls who were well matched for age, but with fewer males and higher IQ. In the ASD group in midline pACC, we found mean 17.7% elevation of glutamate + glutamine (Glx) (p<0.05) and 21.2% (p<0.001) decrement in creatine + phosphocreatine (Cr). We then performed a larger (26 subjects with ASD, 16 controls) follow-up study in samples now matched for age, gender, and IQ using proton magnetic resonance spectroscopic imaging ((1)H MRSI). Higher spatial resolution enabled bilateral pACC acquisition. Significant effects were restricted to right pACC where Glx (9.5%, p<0.05), Cr (6.7%, p<0.05), and N-acetyl-aspartate + N-acetyl-aspartyl-glutamate (10.2%, p<0.01) in the ASD sample were elevated above control. These two independent studies suggest hyperglutamatergia and other neurometabolic abnormalities in pACC in ASD, with possible right-lateralization. The hyperglutamatergic state may reflect an imbalance of excitation over inhibition in the brain as proposed in recent neurodevelopmental models of ASD